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Abstract 

Holas, Howard and March [Phys. Lett. A 310, 451 (2003)] have obtained analytic solutions for 
ground-state properties of a whole family of two-electron spin-compensated harmonically confined 
model atoms whose different members are characterized by a specific interparticle potential energy 
u(ri2)- Here, we make a start on the dynamic generalization of the harmonic external potential, the 
motivation being the serious criticism levelled recently against the foundations of time-dependent 
density- functional theory (e.g. [J. Schirmer and A. Dreuw, Phys. Rev. A 75, 022513 (2007)]). In 
this context, we derive a simplified expression for the time-dependent electron density for arbitrary 
interparticle interaction, which is fully determined by an one-dimensional non-interacting Hamil- 
tonian. Moreover, a closed solution for the momentum space density in the Moshinsky model is 
obtained. 

PACS numbers: 31.15.Ew, 31.25.-v, 31.70.Hq 



1 



I. BACKGROUND AND OUTLINE 



A family of two-electron spin-compensated harmonically confined model atoms has re- 
cently been proposed and studied by Holas, Howard and March |l| (referred to below as 
HHM). Each member of this family is characterized by a specific interparticle potential en- 
ergy w(|ri — r2|). HHM showed that the ground-state spatial wave function \l/(ri,r2) then 
separated in centre-of-mass (c) and relative motion (b) coordinates defined by 

r l+ r 2 K f.s 

c = — - — , b = ri-r a , (1) 

to yield 

*(n, r 2 ) = ^CM(c)^ RM (b). (2) 

The centre-of-mass term ipCM is determined once for all as a Gaussian function. Naturally 
the relative motion part t/'rm entering Eq. (l2"1) , involves the interparticle repulsion u(b), 
but only through an effective one-body potential energy V^(b) given by (here and in the 
following atomic units are used [h, e, m e — 1]) 

V cS (b)= 1 -uj 2 b 2 + u(b), (3) 

which must then be inserted in a one-particle Schrodinger equation to solve for the relative 



motion wave function ipRM- Analytic solutions exist for a number of choices of u(b) in Eq. 
including harmonic [2], Coulombic [3] and inverse square law 0, [|3]. 

As to the motivation for the dynamic generalization of the HHM static family of solu- 
tions in Ref. [lj above, we cite the recent criticism of Schirmer and Dreuw [6] that lies at the 
foundations of time-dependent density-functional theory (TDDFT); a very popular approx- 
imate route for current calculations of electronic excitation energies in atoms and molecules 
fl. Because of such criticisut, « see,„ed the .note hnportant to broadeu —ably the 
class of available exact time-dependent analytical treatments. In that context, we must note 
here for Coulombic interaction in a two-electron model with harmonic confinement [8j], the 
study of D Ami co and Vignale in which the harmonic confinement is generalized to be 
time-dependent, the dimensionality of the problem these authors considered being reduced 
from three in the so-called Hookean atom 3, s| to two. 

Our objective here is to effect a related generalization of the whole family of two-electron 
model atoms considered by HHM to time-dependent theory. We therefore keep the discussion 



as general as possible and turn to specific forms of the interparticle interaction only in 
analytically solvable cases. The resulting expressions could serve as a convenient starting 
point for a simplified numerical treatment, which is however outside the scope of this work. 



II. TIME PROPAGATION OF HARMONICALLY CONFINED TWO-ELECTRON 
MODEL ATOMS 

We wish to solve the Schrodinger equation 



dt 



(4) 



Here we take the model which is a time-dependent generalization of the study of HHM. Thus 
, the Hamiltonian operator assumed has the form 



H = -\ K + V%) + \ul(t) {r\ + rl) + u(|n - r 2 |). 



(5) 



The time-dependent external potential involving cl>q(£) drives the system from its ground 
state \l/o at t=0 to a general time-dependent state \P(ri, r 2 , t) at time t. Following HHM we 
use centre-of-mass vector c and relative motion vector b in Eq. (JT1), to find 



ij^tt(b,c,*) 



*(b,c,t). (6) 



With the help of the product ansatz 



(b, c, t) = V>cm(c, t)^ RM (b, t) 



(7) 



Eq. (jSJ) is without loss of generality readily separated into relative motion and centre-of-mass 
channels according to 



l ^CM(C,t) 



« ^^RM(b,t) 



(8) 



(9) 



-0RM(b,t), 
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with effective masses m cm = 2 and m rm = 1/2. It should be noted, that any (time-dependent) 
separation constant would give rise to phase factors of opposite sign in the RM and CM 
wave functions, which would then cancel in the formation of the total wave function and is 
hence neglected. 



A. Solution of the CM problem in two and three dimensions 

Since the centre-of-mass Hamiltonian does not depend on the interaction potential u(b), 
the CM system be solved once and for all. In two dimensions this has already been accom- 

n 

plished by D'Amico and Vignale [9( for an arbitrary time dependence of uj (t). Their result 
for the CM wave function reads 

^ M > 2D ( C ,t) = J2 C nmXnm(c,t)& m (0), (10) 

n,m 

where the fact that the Hamiltonian does not depend on the centre-of-mass angular variable 
9 was used to separate the wave function into an angular part characterized by the quantum 
number m 

e m (0) = 4= e " me ' (ii) 

V Z7T 



and a radial part 

Xnmict) = A{t)c m e B ^ 2 L™[C{t)c\ (12) 

which involves the generalized Laguerre polynomials L™ and the purely time-dependent 
functions [h]] 



\l - f r, Km<Mi)] 2 
[n + m)\ 



x e -i(2n+m+l)[<P(t)-<P(0)} ^ 



B(t) 
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• /N .d\n\X(t)\ 



(14) 



C(t) = m cm 0(t). (15) 
The complex functions X(t) 

X(t) = \X(t)\e tm with 0(t) > 0, (16) 
are solutions to the equation of motion for the classical harmonic oscillator 

X = -uj 2 (t)X(t), (17) 



which can be solved once an explicit form of uj 2 (t) is chosen. Remarkably, the time prop- 
agation of the harmonically confined quantum system is fully determined by its classical 
analogue. 

As a new result, we compute in the following the CM wave function for the general case 
of three dimensions. The derivation turns out to be straightforward and parallels the 2D 
case with minor modifications. To start with, we write the state ipnim^ with main quantum 
number n, angular momentum I and magnetic quantum number m as a product of a radial 
part Xni and spherical harmonics Yi m : 

^ m (^t) = Xm(c,t)Y lm (0A). (18) 

The more general case in which the system is not in one of its eigenstates at t = can 
be handled easily according to Eq. ffTUj) . since the expansion coefficients do not depend on 
time. Notwithstanding, the ansatz Eq. ( fTBl fully allows for excitations to eigenstates with 
different main quantum number due to the time dependent potential realized by uso(t). 

Proceeding, the radial equation of motion takes the following form: 



.d _ , . 

l-r^Xnl{C,t) 



d 2 Id 1(1 + 1) 1 2 , , 2 

+ 7, o + o^cm^WjC 



Xnl(c,t), (19) 



2m cm dc 2 m cm c dc 2m cm c 2 2 
which we solve using a time-dependent generalization of the well known result for the 
isotropic harmonic oscillator in the ground state 

Xm(c,t) = A(t)^^L^[C(t)(?\. (20) 

Inserting Eq. (1201) into Eq. (lT9"j) and taking advantage of the defining differential equa- 
tion of the associated Laguerre polynomials, we obtain after some algebra the following 
expressions for the functions A, B and C: 

i A + W±VAB-^lAC = (21) 

TTl CTa ^ cm 

ib + —B 2 - \m cm ujl{t) = (22) 

TTl cm Z 

it + —BC + —C 2 = 0. (23) 

Equations fl2TT) to ff23l are (besides different prefactors) completely equivalent to equations 
A4 to A6 in the mentioned work of D'Amico and Vignale, which allows us to write down 



5 



the radial solution in three dimensions immediately: 



Xnl(c,t) 



i 



v/?(n + Z + l)! 



x 



[<j>{t)\ 4 exp 



1 / / / \ .dln|X(t)L 2 
•-m cm (0(t)-z ^^) c 



^(n + f)[0W-0(o)] c / L w/_y mcm ^ (t)c2] . (24) 



B. Time dependent electron density for general interparticle interaction 

Having solved the problem for the CM system we now turn to an evaluation of the time- 
dependent electron density n(r, t) for the special but important case of a system that is in 
its ground state at t — 0. The square modulus of the CM wave function then reduces to a 
simple Gaussian function 



I ,CM,3D/ 
1^000 (c,*)| 



1 



where the time-dependence is fully governed by the characteristic length scale acu.it) of the 
oscillator: 



(25) 



acuit) 



1 



(26) 



m cm 0(t) 

Since the potential u(b) depends on the interparticle distance only, also the relative motion 
wave function ip RM ^ D can once again be split into angular and radial parts according to 



Vw' 3 >,*) = UM)^m(M), 



(27) 



where Cni{b,t) needs in general to be determined by time-propagation in one spatial dimen- 
sion. 

For the mentioned special case the electron density is given by 



n(ri,t) 



'0, 



CM, 3D, 
000 1 



£(ri + r 2 ), t)^' 3D (ri - r 2 , t) 



dr 2 , (28) 



which can be considerably simplified using Eq. fl25|) and fl27|) . After evaluation of the angular 
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integrations and substitution of y = (rj — r 2 )/acM for r 2 , we arrive at 



n(ri,t) 



exp( ^ — — 



x / dy y 2 exp(-^-) ^ooo ' 3D ( a CM(*)y, *) 
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x 



sinh(ri7//a C M(t)) 



(29) 



{r x y / a C yi{t)) 

which constitutes a non-trivial generalization of the HHM result for the static ground state 
density given in Eq. (14) of that publication. 



C. Time- dependent atomic scattering factor 



A quantity which is easily accessible by experiment is the atomic scattering factor given 
by the Fourier transform of the atomic electron density 



/(M) = / n(r,t)e ikr dr 



(30) 



Inserting Eq. into Eq. ( 1301) and taking advantage of the fact that the Jacobian for the 
transformation {ri,^} — > {b, c} is unity, we obtain 



/tat(M) 



^(C^IV^C 



^RM(b,t)|V-rfb 



2/ CM (k,t)/RM(k/2,t), 



(31) 



which shows that the total two-electron scattering factor / to t decouples into a product of 
one-particle scattering factors of the centre-of-mass and relative motion systems. 

While the term /rm remains until the relative motion Schrodinger equation is solved for 
a specific u>o(t) and interparticle potential u(b), the other piece can be evaluated based on 
the results of the last section. In three dimensions a closed solution is however only possible 
for the case of vanishing angular momentum. We therefore step back to the two dimensional 
problem studied by D'Amico and Vignale. For sake of simplicity we treat only the case for 
which the system is in an arbitrary eigenstate at t — 0. The more general superposition of 
Eq. ffTU]) can be solved along the lines of the following derivation and poses no additional 
problems. 
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Using Eq. (T12I) the centre-of-mass structure factor reads 

POO 

/cm(1M) = fcu(k,t) = / J (A;c)|x„ m (c,t)| 2 crfc, (32) 

Jo 

where the angular integration led to the Bessel function of the first kind Jo- After insertion 
of Eq. ({TBI) to (fl5|) and polynomial expansion of the Laguerre functions according to 

s=0 v 7 

we arrive at 

77' 2 

/cm(M) = 7 - r. [m cm <P(t)] m+1 
[n + 171)1 



n n 

X 

s=0 t=0 



EE 



n + m\ /n + 77i\ (— m 



\s+t 



'cms 



n — s J \ n — t J s\t\ 

x / J (A;c)e- mcm * (i)c2 c 2(m+s+t)+1 rfc. (34) 



The solution of the remaining integral is known [Hj and involves the confluent hypergeo- 
metric functions of the Kummer type M: 



POO 

/ e~ aH2 t»- l J v {bt)dt 
Jo 



2a^T(v + 1) 2 " ' ' Aa 2 ' 

V 9fc(// + z/)>0; K(a 2 )>0. (35) 



The CM structure factor then takes the final form 
f C M(k,m c J(t)) = (^yEE (^7) (n-T) 

^ ' s =0 t=0 V 7 V 7 



r(m + s + t + l) „ w _ _ k 2 . , . 
x -i — ^M(m + s + t + 1, 1, r), (36) 



'cms 



which constitutes one of the main results of this work. As indicated, the time-dependence of 
Eq. fl36l) is solely determined by the phase derivative 0, which equals the frequency c^o of the 
confining harmonic potential in the static limit. For a system initially in its ground state, 
/cm reduces to a simple Gaussian which correctly tends, respectively, towards the number 
of electrons (one in this case) as k approaches zero, and zero as k grows to infinity. 

Having obtained a closed solution for the centre-of-mass structure factor we now briefly 
discuss the relative motion term. In general it will be necessary to evaluate this part by 
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numerical methods, which is beyond the scope of the present article. The so-called Moshin- 
sky atom [2| characterized by the harmonic interparticle potential u(b) = — \K>b 2 provides 
a useful exception and can be treated analytically. To this end we note that in this case 
the relative motion Hamiltonian is equivalent to a centre-of-mass one with effective mass 
m cm = 1/2 and force constant ^(t) — /C/m cm . Solving 

X(t) = \X(t)\e i ^ t) with X = - (ujlit) - — ^ X, (37) 

V m cm J 

allows one to obtain the total structure factor of the Mochinsky atom from 

f t K ot (k,t) = 2f C M(k,m cm ^t))f cu (k/2,m cm 4>(t))- (38) 

III. RELATION TO TDDFT 

As mentioned in the Introduction, this work was motivated by the serious criticism of 
Schirmer and Dreuw [(J of the 1984 theorems of Runge and Gross [ijj]. In Ref. 6] the proof 



m 



12) wa 8 not only changed bnt seemingly refuted. To be speeifle, the anthers elahned 
that the variational derivation of the time- dependent Kohn-Sham equations in Ref. [12J is 
invalid due to an ill-defined action functional presented there. A nonvariational formulation 
would also run into problems, since in this case the Kohn-Sham system would allow one to 
reproduce but not to predict the exact electron density. 

The results of this work can in principle be used to investigate this issue in an actual 
numerical calculation. For the two-electron spin-compensated system at hand, it is possible 
to construct the exact exchange-correlation potential from the known electron density in 



real (Eq. [291) or momentum space (Eq. I3B1 . This is in line with van Leeuwen's proof 13] of 



mapping from densities to potentials in TDDFT with general two-particle interactions. Such 
a construction was for example already performed by Lein and Kiimmel or D Amico and 
Vignale jij . Using this potential in the numerical solution of the time-dependent Kohn-Sham 
orbitals would then open the opportunity to compare the propagated and exact electron 
density. 



IV. SUMMARY AND POSSIBLE FUTURE DIRECTIONS 

As set out in the Introduction, HHM [1] proposed and solved the problem of a two- 
electron spin-compensated family of harmonically confined atoms with a general interparticle 



repulsive potential u{r\2)- 

Using the available time- dependent study of D'Amico and Vignale jjj, but now in two 
dimensions, we display here the time-dependent density n(r, t) for a system initially in the 
ground state in three dimensions. While the general result is somewhat formal, for the 
so-called Moshinsky atom characterized by the choice u(b) = — |/C6 2 with K measuring 
the strength of the interparticle repulsion, we derive the momentum space electron density 
f K (k,t) functionally from the D'Amico and Vignale form for the independent two-electron 
case with /C set equal to zero. 

For the future, progress may come by taking a specific choice of the time-dependence oj^it) 
of the harmonic confinement potential. Singling out a specific Fourier component through 
the choice expitut+exp — iujt is restrictive to periodic time-dependence but may allow, in the 
future a more elegant discussion of the time-dependent particle density n(r,t) at the heart 
of all current theories of TDDFT. But, of course, more important would be to pass from 
a model 'He-like' family to a four-electron interacting 'Be-like' system, where an elegant 



Dirac idempotent density matrix already exists in 



or Kohn-Sham) limit (see, for example Refs. 15l. Il6|). 



;he independent particle (Hartree-Fock 
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